High Reynolds number magnetohydrodynamic turbulence using a Lagrangian model 
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With the help of a model of magnetohydrodynamic (MHD) turbulence tested previously, we 
explore high Reynolds number regimes up to equivalent resolutions of 6000 3 grid points in the 
absence of forcing and with no imposed uniform magnetic field. For the given initial condition chosen 
here, with equal kinetic and magnetic energy, the flow ends up being dominated by the magnetic 
field, and the dynamics leads to an isotropic Iroshnikov-Kraichnan energy spectrum. However, 
the locally anisotropic magnetic field fluctuations perpendicular to the local mean field follow a 
Kolmogorov law. We find that the ratio of the eddy turnover time to the Alfven time increases with 
wavenumber, contrary to the so-called critical balance hypothesis. Residual energy and helicity 
spectra are also considered; the role played by the conservation of magnetic helicity is studied, and 
scaling laws are found for the magnetic helicity and residual helicity spectra. We put these results 
in the context of the dynamics of a globally isotropic MHD flow which is locally anisotropic because 
of the influence of the strong large-scale magnetic field, leading to a partial equilibration between 
kinetic and magnetic modes for the energy and the helicity. 

PACS numbers: 47.65.-d, 47.27.Jv, 94.05.Lk, 95.30.Qd 



> 
00 

in 
in 

o 



I. INTRODUCTION 

Magnetic fields are known to be dynamically impor- 
tant in a variety of flows, for example in the liquid core 
of the Earth, in stars and galaxies, and they also play 
a role at the atomic level, e.g., in Bose-Einstein conden- 
sates [l|. Their origin in astrophysics and geophysics is 
still ill-understood, and their dynamics in the fully non- 
linear regime once the magnetic field has grown, possibly 
to quasi-equipartition with the velocity, is ill-understood 
as well. Numerous studies have shown, though, that mag- 
netohydrodynamic (MHD) turbulence differs in several 
ways from classical hydrodynamic turbulence; for exam- 
ple, it is more intermittent than fluids [l-Q, more so in 
particular in two space dimensions, and this intermit- 
tency may be either time-dependent or dynamically con- 
strained as observations of solar active regions show [f| . 
Furthermore, numerical and theoretical arguments indi- 
cate that MHD turbulence is less spectrally local than 
hydrodynamic turbulence 0, 0] (see also P-fToj). and 
that MHD may lack universality in its scaling properties 

EMU, 

including in the absence of an imposed uniform 
magnetic field [15| . Finally, the role of the correlation be- 
tween the velocity and the magnetic field has been known 
for a long time to be important, and this was particularly 
emphasized when deriving an exact law involving linear 
scaling with distance of third-order structure functions 
[lH (see [TtJ for the von Karman relations for MHD in 
terms of correlation functions) . 

Part of these differences may stem from the fact that, 
even in the incompressible MHD case, the system sup- 



ports waves due to the effect of a large-scale magnetic 
field (hence the non- locality), whereas there are no such 
waves in incompressible hydrodynamics in the absence 
of external agents such as gravity or rotation. Of course, 
at small scales, plasma kinetic effects will come into play 
(e.g., in space plasmas), rendering the analysis more diffi- 
cult. But even when one ignores such small-scale effects 
and concentrates on the dynamics of the large scales, 
one is still left with a variety of observations and of phc- 
nomenological approaches to the problem. 

A wealth of new space data has come with the Cluster 
mission (which uses four identical spacecrafts to study 
the Earth's magnetosphere) , and future projects such as 
the Magnetospheric Multiscale (MMS) mission will allow 
for more in-depth studies of MHD turbulence, albeit in a 
complex setting involving compressibility, boundary con- 
ditions, and geometrical effects as well as the aforemen- 
tioned plasma phenomena. MHD turbulence can also be 
studied in the laboratory, although achieving high mag- 
netic Reynolds numbers in liquid metals will remain a 
challenge for some time to come. 

Of particular interest is how energy is distributed 
across scales in MHD turbulence. For the moment there 
is no clear-cut answer, with some indications that uni- 
versality may be breaking [l5| . but with the parame- 
ter (s) distinguishing between possible different classes of 
solutions not well-known. Direct numerical simulations 
(DNS) played an important role in these findings, as they 
allow for a rather controlled exploration of solutions, with 
well-defined parameters, but with the disadvantage that, 
in three space dimensions, the resolutions that are at- 
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tainable with present-day computers are insufficient, by 
far, if one wants to model realistic geophysical and as- 
trophysical flows with huge Reynolds numbers. Also, the 
computational cost of a detailed study of parameter space 
at sufficiently high resolution is currently out of reach. 

One way out of this conundrum, in part, is to resort 
to modeling. There are numerous methods that have 
been devised over the years (see, e.g., recent reviews for 
fluids [Hi and for MHD \v^). Among them, the La- 
grangian averaged MHD model (LAMHD hereafter) de- 
veloped in [2(| (see also [2l|) seems promising in that 
it allows to perform a quasi-DNS, in the sense that the 
Reynolds number is known and that the inertial range 
is extended, compared to a DNS performed on the same 
grid without the model, thanks to a different formulation 
of the equations that preserves the invariants, albeit in a 
different norm (sec below). 

In this context, we investigate properties of decaying 
MHD turbulence in the absence of an imposed uniform 
magnetic field using the LAMHD model. We thoroughly 
tested the LAMHD model previously, in two dimensions 
0, HH, [HI an d in three dimensions (23J. Here we ex- 
tend our computations to an equivalent grid resolution 
of the order of 6000 3 points, with Taylor Reynolds num- 
ber at the maximum of dissipation of w 2900. The next 
section is devoted to definitions and the methodology, 
Sec. Mil gives details of the simulation, Sec. IIVI discusses 
the global dynamics, and Sec. [V] discusses the role played 
by Alfvenic exchanges. Finally, Sec. lVIl is our conclusion. 

II. EQUATIONS 

The MHD equations for an incompressible fluid with 
v and b respectively the velocity and magnetic fields in 
dimensionlcss Alfvenic units are: 

<9 t v + wxv=jxb - Vp + fV 2 v , 
d t a = v x b - + r/V 2 a , 

V-v = , V-a=0, (1) 

with b = V x a, a being the magnetic vector potential 
in the Coulomb gauge. The potential (f> and the pressure 
divided by the constant (unit) density p are obtained self- 
consistcntly to ensure respectively the Coulomb gauge 
and the incompressibility of the velocity field. The 
kinematic viscosity is v and the magnetic diffusivity 77; 
uj = V x v is the vorticity, and j = V x b is the cur- 
rent density. The magnetic Prandtl number Pm = v/f) 
is taken in what follows equal to unity. 

In the absence of dissipation (y = , r\ = 0), the total 
energy 

E T = E V + E b = -(||v|| 2 + ||b|| 2 ) = 

= \ J (|v| 2 + |b| 2 )d 3 x, (2) 

is conserved (E v and E^ are respectively the kinetic and 
magnetic energy, and || • || 2 denotes £ 2 norms). Two 



other quadratic quantities are preserved by the nonlin- 
ear terms in three dimensions (3D): the cross-correlation 
He = J v • b d 3 x, and the magnetic hclicity Hi, = 
jab d 3 x. 

The Lagrangian averaged model introduces a priori two 
filtering lengths, a„^, for the velocity and the magnetic 
fields in such a way that the conservative structure of the 
equations is preserved (20| . Written in terms of the rough 
fields v, b, and a, and of the filtered fields u, B and A, 
the LAMHD equations read: 

a ( v + uxu = jxB - VII + i/V 2 v 
d t A = uxS - V</>' + ??V 2 a 
V-v = V-u = V-a=V-^ = 0, (3) 

where II is a modified pressure, <fi' a modified poten- 
tial, and with filtering being accomplished through nor- 
malized convolution filters L Vj b chosen to be the in- 
verse of a Hclmholtz operator, namely L v ^ = H^ 1 = 
(I — a 2 h V 2 ) -1 . The filtering lengths are related to 
the ratio of the dissipation scales of the LAMHD flow 
and the modeled MHD flow [24|; since we shall investi- 
gate flows with unit magnetic Prandtl number, we take 
a common filter length for both the velocity and mag- 
netic field, a Vt b = a. That is, u = (1 — a 2 V 2 ) _1 v and 
B = (I — a 2 V 2 ) _1 b. The ideal quadratic invariants in 
terms of £ 2 norms for MHD arc transformed for LAMHD, 
being expressed now in terms of H" norms; this leads to 
formulations that contain both the rough fields v and b 
and the filtered fields u and B, as can be seen for example 
in the expression of the total energy: 

K = ^(l|u||f + l|B||f) 

= 2/ [( u -" 2 V 2 u) -u+ (B-a 2 V 2 B) -S] d 3 x 
= - J [v u + b • B] d 3 x. (4) 

Similarly, the magnetic helicity H" = J a • B d 3 x in its 
Hi norm is preserved by the LAMHD equations. A form 
of cross-correlation is also preserved, and the LAMHD 
equations also have the equivalent of an Alfven flux- 
conservation theorem (based on these equivalences, and 
for the sake of clarity, the superscript a will be suppressed 
in the sections presenting numerical results unless strictly 
required) . 

The Reynolds number in the following is defined as 

Re = U rms L /u , (5) 

with U rms and L the r.m.s. velocity and the integral 
scale respectively. The Taylor Reynolds number R\ = 
U rms X/v is defined using either the kinetic or magnetic 
Taylor scale: 

Kb = ^y/Kb/Kb . ( fi ) 
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where fi" = and = = ||j||2 are respectively 

the enstrophy and square current (or magnetic enstro- 
phy) in LAMHD, proportional respectively to the kinetic 
and magnetic dissipation in the model. Taylor wavenum- 
bers are defined as k\ v b = 2n/X v ^- 

The numerical implementation of the LAMHD equa- 
tions is performed using a pseudo-spectral code for which 
inversion of the Helmholtz operator is straightforward. 
The code is parallelized using MPI (Message Passing 
Interface) and has been tested in a variety of condi- 
tions (see, e.g., [HJ]); a hybrid version of the code, using 
OpcnMP as well, allows efficient parallelization for higher 
resolutions using a larger number of processors [26j . The 
computational box has a size of (27r) 3 , and wavenumbers 
vary from k m i n = 1 to k ma x = N/3 using a standard 2 /3 
de-aliasing rule, N being the number of grid points per 
direction. 



III. SIMULATION SET-UP 

The dynamics of MHD turbulence is quite complex and 
there may be different regimes, as characterized for exam- 
ple by their energy spectra, either isotropic or anisotropic 
depending on the presence of an imposed uniform mag- 
netic field bo. With bo = and no forcing, a high reso- 
lution DNS using 1536 3 grid points [27| showed evidence 
of isotropic Iroshnikov-Kraichnan (hereafter, IK) energy 
scaling [28||. More recently, it was shown [l5[ that en- 
ergy spectra may differ measurably even when global 
quantities (Et, E b /E v , H b , and He) are the same at 
the initial time: IK-like, weak turbulence like (hereafter, 
WT) [29|], or Kolmogorov-like spectra (hereafter, K41) 
[30j were obtained, the discriminating factor being the 
ratio of Alfvenic time ta to the eddy turn-over time tnl ■ 
These times are defined respectively as: 



ta = e/b 



TNL = £/v£ , 



ve being the characteristic velocity at scale £; in the ab- 
sence of an imposed mean field, bo in the above expres- 
sion is taken as the amplitude of the magnetic field at 
the gravest mode. 

In this paper, we examine one of these cases further 
by performing computations at higher Reynolds number 
using the LAMHD model, resulting in an equivalent grid 
of roughly 6000 3 points and a Taylor Reynolds number 
of w 2900. The case of interest is the one that gives IK- 
like energy scaling, and we therefore consider only one 
initial condition, similar to the one studied in Ref. (27j 
with 1536 3 spatial resolution. By "equivalent grid," we 
mean the following: a comparison with the 1536 3 DNS 
with r/ = is = 2 x 10~ 4 was carried solving the LAMHD 
equations on grids of 256 3 , 384 3 , and 512 3 points to vali- 
date the results of the model [Til HH . The runs presented 
here, with 1024 3 grid points and viscosity and magnetic 
diffusivity ?] = z/ = 5xl0~ 5 , would require a « 6000 3 
resolution in a DNS. 



The filter length a is chosen to be ^/k max , i.e., the 
filter wavenumber is fc maa ,/2, and the grid resolution is 
Ax = 27rfc~i, r /3 (see, e.g., [23| for details on how the 
filter length is chosen for a given resolution). The initial 
conditions are close but not identical to those used in 
[27J ; they consist of a superposition of Arn'old-Beltrami- 
Childrcss (ABC) flows distributed in the first four shells 
(k G [1,4]), with super-imposed noise with a wide spec- 
trum as in [27| . However, the noise has a steeper energy 
distribution, oc fc -6 , instead of oc fc~ 3 . The choice of a 



shallower distribution of random noise in [27j was guided 
by the desire to reach the maximum of dissipation in a 
short computational time; the computational constraint 
not being so strong in the LAMHD model, we resort to 
a steeper initial random distribution so as to let turbu- 
lent spectra evolve in a less constrained manner, allowing 
us to examine in this paper as well whether qualitative 
differences in behavior take place or not under such a 
change. 

The steeper distribution of noise here also results in 
slightly different values for the relative initial helicities 
(when compared with the run in (27j). defined respec- 
tively as 



p v = H v /^/E v fl v 



Pb = H b /y/E a E, 



Here, E a is the £2 norm of the vector potential, and 
H v = J u • u) d 3 x is the kinetic helicity, which measures 
deviations from mirror symmetry in the hydrodynamic 
flow. Because of Schwarz inequality, we have \p v ,b\ < 1- 

The resulting flow at t = has E v = E b = 1, p v = 
0.78, and p b = 0.94, compared to E v = E b = 1, p v = 
0.49, and p b = 0.54 for the DNS run in [27|, leading to 
stronger nonlinearitics and faster initial evolution in the 
previous work. Note that at t = we can use either £2 
or Hi norms for the LAMHD energy and other global 
quantities, since they differ only after the 4th digit. 



IV. GLOBAL DYNAMICS 

A. Time evolution 

As expected from the larger values of relative helicities, 
development of a fully developed turbulent regime takes 
a longer time in the present run compared to the DNS 
run on a grid of 1536 3 points. Indeed, the maximum 
of total dissipation for the LAMHD flow is reached for 
U w 6.5 (as opposed to w 3.8 for the 1536 3 DNS run) 
as can be seen in Fig. Q] (inset). The total dissipation 
remains approximately constant between the maximum 
and the final time of the computation, t pa 11, providing 
us with a quasi-steady state at small scales for several 
turn-over times. Except for this delay in the onset of 
the formation of small scale gradients, both flows evolve 
in similar manners. Figure 1 also shows the temporal 
evolution of the total energy, and of the ratio of magnetic 
to kinetic energy. The magnetic energy grows rapidly 
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FIG. 1: Ratio of total magnetic energy to kinetic energy 
E" as a function of time, and (inset) total energy, E% = 
E™ + E£ (solid line), and total dissipation, e a = vQ.^ + r)Q, b 
(dash-dotted line) versus time. Note that the superscript a 
is suppressed in the labels, and will be suppressed in later 
figures for clarity. 



at the expense of its kinetic counterpart, and the ratio 
between these two quantities then fluctuates around a 
value close to « 3.5 for times between » 3 and « 10, 
with a tendency for further increase toward the end of 
the run. 
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B. Energy scaling 

During the time of approximately constant dissipation 
rate, the energy spectrum presents a clear IK law, as can 
be seen in Fig.[5Ja), which displays the total energy spec- 
trum compensated by fc 3 / 2 . In this figure, and in the sub- 
sequent ones, the vertical solid line represents the classi- 
cal Taylor scale based on the velocity field, the vertical 
dash line the magnetic Taylor scale, whereas the verti- 
cal dotted line indicates the filter wavenumber 2n/a (be- 
yond which the model is expected to differ from MHD). 
Note that no bottleneck is observed, i.e., there is no en- 
hancement of the spectrum at the onset of the dissipa- 
tive range, as already observed in several MHD simula- 
tions. The IK law seems a good approximation to the 
dynamics in the range of wavenumbers from fc « 4 to 
k « 27r/Af, « 28, i.e., all the way to the magnetic Taylor 
scale. 

To verify this, Fig. [2jb) also shows the total, kinetic, 
and magnetic energy spectra compensated by fc 19 / 12 . 
The spectral index 19/12 corresponds to the mean value 
between the IK and K41 spectral indices, i.e., (3/2 + 
5/3) /2. As a result, a positive slope in the compensated 
plot indicates the spectrum is closer to IK, while a nega- 
tive slope indicates the spectrum is closer to K41 scaling. 
The total and magnetic energy spectra monotonously in- 
crease up to the magnetic Taylor wavenumber, and only 



FIG. 2: (a): Er(k) averaged over time witht £ [6.5, 11], com- 
pensated by fc 3//2 (solid line, left scale) and ratio of magnetic 
to kinetic energy as a function of k ( dash-dotted line, r ight 
scale). Also shown is -Er(fc) from a lower Re 1536 3 DNS \2m 
(dashed line) demonstrating the efficacy of the model. Here 
and in all subsequent figures, the vertical solid, dashed, and 
dotted lines denote respectively the wavenumbers correspond- 
ing to the kinetic Taylor scale, the magnetic Taylor scale, and 
the filtering scale a (beyond which the model is expected to 
differ from MHD with n = v = 5 x W~ 5 ). (b): Total (solid 
line), magnetic (dashed line) and kinetic (dotted line) energy 
spectra compensated by k 19 ^ 12 = k^ 3 ^ 2+3 ^ 3 ^ 2 (see text for de- 
tails). The IK slope (—3/2) is favored against the K41 (—5/3 ) 
slope until the magnetic Taylor scale. 

the kinetic energy spectrum flattens and has a small 
negative slope between the kinetic and magnetic Taylor 
wavenumbers. 

Figure[2fa) also gives (with scale on the right) the ratio 
E"(k)/E"(k); it appears remarkably flat starting at the 
beginning of the IK incrtial range, and it also displays 
a slight enhancement peaking around fc = 80 (associated 
with the faster decay of the kinetic energy spectrum after 
the kinetic Taylor scale). The approximate constancy, 
with Eb(k) m 1.4 E v (k), denotes a partial Alfvcnization 
of the flow in most of the inertial range, but with here 
at k = k m i n a magnetic field which can be evaluated as 
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FIG. 3: Residual energy spectrum En(k) = E b (k) — E v (k) 
averaged in the same time interval as in Fig. [2] and compen- 
sated by fc 7 / 3 . The dashed line indicates a slope of —2 up 
to the Taylor magnetic scale, while the solid horizontal line 
indicates a slope of —7/3; the best fit from 4 < k < 28 gives 
-1.4 ±0.1. 
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FIG. 4: Ratio of the eddy turn-over time tnl to the Alfven 
time ta as a function of wavenumber averaged in the same 
interval of time as in Fig. [2] (solid line); the power law fit 
(dot-dash line) has an exponent of 0.3 with an error of ±0.1. 



being roughly 3 times larger than the large scale velocity 
field, and at least 10 times larger than the r.m.s. value of 
the turbulent fluctuations. 

Note that, at scales smaller than the filter a, the bal- 
ance between the velocity and the magnetic field for the 
model is inverted, with a slight excess of kinetic energy 
(e.g., note E"(k)/E"(k) < 1 in this range). This may be 
attributed to an effective hyper-resistivity in LAMHD for 
the magnetic field in the sub-filter scales; at these scales, 
the model differs from MHD. Indeed, in the DNS re- 
ported in [l5| , this decrease of Et /E v was only seen close 
to the dissipation frequency and the magnetic to kinetic 
ratio remained above unity at all scales. 

The residual energy spectrum 

E R (k)=EZ(k)~E:(k) : (7) 

is shown in Fig. [3] A heuristic argument indicates that 
for K41 scaling, the residual energy is expected to vary as 
fc -7 / 3 , while for IK scaling the residual energy should go 
as fc -2 (32|. In the inertial range of our LAMHD run, the 
slope is closer to (although shallower than) —2 (a result 
also seen in Ref. j33[.) This point is further discussed in 
Sec. V. 

We show in Fig.U]the ratio t tvl(^)/'7~a(^<) as a function 
of wavenumber, where the eddy turnover time and the 
Alfven time are defined as 

r NL (k) = — - 1 (8) 

TA{k) = \jk*Ez(k = iy (9) 

respectively. 



Since Eg(k)/E£(k) ~ constant (see Fig.^a)), this im- 
plies that the associated characteristic time ratio ta/tml 
must be a function of wavenumber. Indeed, the ratio is 
clearly not constant as the wavenumber increases, and 
seems to follow in the inertial range a k 3±01 power 
law (Fig. g]). If we assume E v (k) ~ E T (k) ~ fc" 3 / 2 , 
this leads to tjvl ~ fc~ 3 / 4 and to the ratio shown in 
Fig. 0] varying as fc 1 / 4 . When we assume K41 scaling 
with E v (k) ~ fc~ 5 / 3 , the ratio of times varies as fc 1 / 3 . 
Finally, for WT we would have fc 1 / 2 . The first two so- 
lutions are compatible with the data, in particular con- 
sidering we are not taking into account that (i) intcrmit- 
tency corrections to the spectral indices can alter these 
indices, and (ii) the kinetic and total energy spectra may 
not follow the same power law. Indeed, one should re- 
member that there is no constraint on the behavior of 
E v fi{k) separately, since the invariant that cascades di- 
rectly is the total energy. In fact, the stronger constraint 
on the spectral behavior of MHD fields comes from the 
exact law in terms of the energy and energy fluxes of 
Elsasser variables, or in terms of the total energy and 
cross helicity and their fluxes (as well as for magnetic he- 
licity). These laws imply complex correlations between 
the velocity and the magnetic field [ljj], but do not im- 
ply that power laws must be followed by each individual 
field. Finally, note that this result is not in agreement 
with the critical balance hypothesis [3l[ which postulates 
TNLik) = TA{k) ~ fc -1 ; in fact, solving this relationship 
for the kinetic energy spectrum E v (k) this hypothesis 
would lead to E v (k) ~ fc -1 in the isotropic case, imply- 
ing a logarithmic divergence which is unlikely (the case 
with an external uniform magnetic field may be differ- 
ent). 
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FIG. 5: Spectrum of magnetic helicity averaged over the same 
interval of time as in Fig. [2] (solid line); the dot-dashed line 
indicates a power-law fit with exponent —3.7 and an error of 
±0.1. 
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FIG. 6: Residual helicity H R (k) = k 2 H b (k) - H v (k). The 
dashed line indicates a power law ~ k~ , while the dash- 
dotted line indicates ~ k~ 2 ' 2 . 



C. Helicity scaling 

Examining now the magnetic helicity spectrum, we ob- 
serve that it shows a single scaling in the energy contain- 
ing range, in the inertial range, and beyond the Taylor 
scales (in contrast to the spectra discussed before). It 
follows rather closely a fc- 3 - 7 ^ - 1 law, as can be seen in 
Fig. [5j Such a steep spectrum for Hb has been observed 
before 0, HE] . The slope does not correspond to the phe- 
nomenological analysis done in [32l ] under the hypothesis 
of an inverse cascade of magnetic helicity. Indeed, in 
this simulation all the helicity is concentrated initially in 
k 6 [1,4], and therefore there is no separation between 
the helicity containing scale and the size of the box for 
an inverse cascade of magnetic helicity to develop. The 
power law in Hb(k) seems therefore to be associated with 
a transfer of the helicity towards smaller scales by the di- 
rect energy cascade. As will be discussed in the next sec- 
tion, this transfer is required to satisfy a balance imposed 
by Alfvcnization and magnetic field induction. Note that 



34[, such a steep spectrum for Hb(k) is interpreted as 



a Kolmogorov spectrum for both the current helicity and 
the kinetic helicity, through sweeping by the large scales. 

Similarly to the case of energy, a residual helicity spec- 
trum can be defined as 



H R (k) = k 2 H b (k)-H v (k). 



(10) 



The factor k 2 multiplying Hb(k) gives the current helicity 
spectrum, which is dimensionally equivalent to H v {k). 
This functional is known to play an important role in 

Fig. [5] shows 



32] 



the nonlinear stage of the dynamo 
Hn(k) in the simulation. Two ranges become apparent, 
with a knee in the spectrum around k Ri 9. For larger 
wavenumbers, the residual helicity seems to follow a ~ 
k- 2 - 2 law. We observe that H b (k) > Vfc, and k 2 H b (k) - 
H v (k) is positive until k « 100. 



D. Local anisotropy 

We further examine in Fig. [7] the local scaling of the 
magnetic field when expressed in terms of its second- 
order longitudinal structure function, separated into per- 
pendicular and parallel spatial increments, 



s b 2 (h 



[b(x+l ± ,||)-b(x)] 




(11) 



Here, _!_ and || denote space-varying directions with re- 
spect to the direction of a local mean magnetic field 
bo.iocj defined as an average of b in a box of edge length 
7r (the flow integral scale is Lq w 2.9) around each point 
x for which the data is sought. Therefore, these struc- 
ture functions do not measure global anisotropy (the flow 
is globally isotropic) or global scaling, but rather local 
anisotropy and scaling with respect to a local mean mag- 
netic field. 

The points at which these correlation functions are 
computed are chosen at random, and the figure is the 
result of an analysis performed for 2.5 x 10 7 points. Note- 
worthy is the variation with scale of the ratio of the 
perpendicular to the parallel component, with a cross- 
over near the magnetic Taylor scale. At smaller scales, 
S^J-) > £2(^)1 an d the ratio increases as the scale de- 
creases until reaching the filter scale. This is in sharp 
contrast with hydrodynamical turbulence, including in 
the presence of an imposed rotation: isotropy is recov- 
ered at small scale in such flows because the ratio of the 
inertial wave time to the eddy turn-over time gets larger 
as the wavenumber grows, or in other words the influence 
of anisotropy is felt at large scale but becomes secondary 
at small scales. It is also noteworthy that the local per- 
pendicular magnetic field seems to display a scaling closer 
to K41 turbulence, while the parallel increments, having 
SjjGll) ~ 'ih display a scaling close to ~ k7 2 . 
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FIG. 7: Second-order longitudinal structure function (for 
t = 7) of the magnetic field, compensated by a Kolmogorov 
law Z 2 / 3 , for perpendicular (solid line) and parallel (dotted 
line) increments with respect to the local mean magnetic field. 
The structure function for parallel increments (dotted line) 
displays, at scales larger than the Taylor scales, a scaling ~ In. 



V. ALFVENIC EXCHANGES 

Is there a dynamical role for Alfvenic exchanges be- 
yond a tendency toward equipartition? Even if the spec- 
tral behavior of the total energy can be obtained from 
first principles and phenomcnological considerations as 
to what is the relevant timcscalc for energy transfer (the 
eddy turn-over time, the Alfvcn time, or their combina- 
tion in a transfer time r tr = t^ l /ta that embodies the 
slowin g d own of nonlinear coupling in the presence of 
waves [28|), the behavior of individual spectra is not so 
straightforwardly determined. For the kinetic and mag- 
netic energy spectra, there are indications of differing 
spectral behavior both in the solar wind [HI and in nu- 
merical simulations 



37j. For the magnetic helicity, flows 



forced at intermediate scales have been observed to un- 
dergo an inverse cascade to large scales including in the 
supersonic case (see [HI and references therein), but in 
the decaying case and with initial conditions at large 
scale, the behavior of H b {k) does not appear to be uni- 
versal [35| . In this section, we discuss the scaling laws ob- 



served in our simulations in terms of Alfvenic exchanges 
and magnetic induction. 

We start with the residual energy En(k), Fig.j3l The 
way to derive En(k) ~ fc~ 2 is to suppose that [33]: (i) 
the residual energy is small compared to the total energy, 
in a ratio (t^/ttvl) 2 , i.e., 



E R (k) ~ (T A /T NL ) 2 E T (k) 



(12) 



and (ii) that the kinetic energy spectrum appearing in 
t~nl (see Eq. [5]) behaves in the inertial range as the total 
energy. Assuming further that Et follows an IK k~ 3 ^ 2 
law leads to Eji(k) ~ fc~ 2 , but the second hypothesis 
may not be fulfilled in all cases [HI, [13] and indeed we 



find a shallower residual spectrum, Er ~ fc -3 / 2 . 

The understanding of the origin of the discrepancies in 
the scaling of magnetic helicity, Fig. [5J is partial at best. 
We focus here on the spectral scaling in the direct cas- 
cade energy range presented in the previous section, as 
opposed to the inverse cascade of magnetic helicity which 
is not accessible in our simulations. It is advocated in [35| 
that the scaling can be attributed to a combination of fac- 
tors involving the competing nonlincarities in MHD tur- 
bulence, modeled using the two-point closure developed 
in (32l ] . A phenomenological way to recover the balance 
advocated in [35| is to assume partial Alfvenization of 
the flow, leading to a quasi-equipartition of the velocity 
and magnetic field, both for the symmetric part (ener- 
gies) and the anti-symmetric part (helicities) of their re- 
spective correlation tensors. In other words, one could 
expect to have, in the order of an Alfven time (based on 
the large-scale magnetic field), the following balance 



E v (k) ~ E b (k) , H v (k) ~ k 2 H b {k) 



(13) 



Note however that equipartition is not achieved exactly, 
as observed in the solar wind. One can then postulate 
that this incomplete equipartition be equally occurring 
for the energy and for the helicity. That is, that the 
ratio 



E v (k)/E b (k) 
H v (k)/[k 2 H b (k)} 



(14) 



be close to a constant in the inertial range. This condition 
is approximately fulfilled in our simulation as well as in 
the DNS in [13,|33. 

The relationship in Eq. (|14[) can be derived in a differ- 
ent manner, that takes into account the dynamics of the 
magnetic field and further clarifies how the system can 
evolve towards such balance. To do this, we take a mean- 
field approach to induction processes in the equation for 
the magnetic field, associated in the mechanically forced 
case with the growth of a seed magnetic field by dynamo 
action. By splitting the magnetic field into b = bo + b', 
where bo is a large-scale magnetic field that grows un- 
der the influence of the small-scale kinetic helicity, and 
b' are the turbulent magnetic fluctuations, one obtains 
under the mean-field approximation: 



db 

at 



a'V x b . 



(15) 



Here a' oc — r corr (H v ) ss , where (.) indicates averaging 
over the small-scale turbulent helical field, and T corr is 
a velocity field correlation time [3!| (see [ioj] for a re- 
view). Taking the dot product of the above equation 
with the large-scale magnetic potential a and integrat- 
ing over space, one gets the following evolution equation 
for the large-scale magnetic helicity: 



dt 



a • h d 3 x = 2a' / |b | 2 d 3 x 



(16) 



Note that dissipation has been neglected in these 
equations, and that the large-scale magnetic energy 
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(J | bo \ 2 d 3 x) appears in the equation upon integration un- 
der suitable boundary conditions. 

As a result of Eq. (fT()|) . the helical velocity field injects 
helicity into the large-scale magnetic field. Since total 
magnetic helicity is conserved in the ideal case, the op- 
posite amount of magnetic helicity must be created at 
small scales, i.e., the following dynamical balance must 
hold independently of the scale (with a change in sign for 
large and small scales): 

dHh 

—— ~ ±T corr H v Eb . (17) 
at 

The dynamical balance can be understood using the con- 
ceptual "stretch, twist, and fold" (STF) dynamo [Zlj . 
Each time a closed magnetic flux tube is twisted by the 
helical velocity field, smaller-scale magnetic field lines are 
twisted in the opposite direction, thus pushing towards 
small scales some magnetic helicity. This process re- 
moves magnetic helicity from the large scales and allows 
the magnetic field to "disentangle" through reconnection 
events, destroying in that way magnetic helicity fiol l42j . 

The correlation time r corr can be taken to be the 
eddy turn-over time, r corr ~ t^l, since the induction 
effect given by Eq. ([17]) is associated with deforma- 
tion of magnetic field lines by turbulent eddies. From 
Eq. (jTTf , assuming again that the temporal evolution fol- 
lows the eddy turn-over time tjvx, i-e., by approximating 
d/dt ~ 1/tjvl, we find 

H b oc H v E b /(k 2 E v ) (18) 

from which Eq. (frjj follows. These phenomenological 
considerations show that the regeneration of the mag- 
netic field by dynamo action is in fact compatible with 
an Alfvenic balance, and that these two processes may 
be occurring on similar time scales in MHD turbulence. 
In a decaying flow such as here, there may not be suf- 
ficient time available to see such Alfvenic and dynamo 
exchanges take place repeatedly, but it would be of in- 
terest to examine the helicity exchanges in a forced flow; 
this may provide a geometrical origin of Alfvenization, 
through the dynamic alignment of relevant vectors, as is 
known to occur on a nonlinear time scale (43j . 

The constraint given by Eq. (|18l) is compatible with 
the spectrum of magnetic helicity observed in Fig. [5] In 
hydrodynamic turbulence, k~ 1 H v jE v ~ 1/k [441 ) . and 
under this assumption Hb ~ k~ 2 Eb (in this simulation, 
the ratio k~ 1 H v / E v decays slightly faster with k). For 
Eb ~ Et ~ fc~ 3 / 2 or steeper, then Hb ~ fc~ 7 / 2 or steeper. 
The fc- 3 - 7±01 law in Fig. seems therefore consistent 
with a transfer of the helicity towards smaller scales by 
the direct energy cascade, to satisfy the balance imposed 
by Alfvenization or induction. 

We finally consider the spectrum of residual helicity 
from the point of view of Alfvenization. To derive its 
scaling, we assume that, similarly to the residual energy 
case, we have: 

H R (k) ~ (T A / TNL ) 2 H c (k) . (19) 



The cross-correlation appears here for at least two rea- 
sons: On the one hand, helicity is a pseudo-scalar and 
thus the r.h.s. of Eq. (jT9"]) must involve a pseudo-scalar 
as well. On the other hand, it is known that total en- 
ergy and cross-correlation play inter-linked roles in the 
dynamics of MHD flows. This is particularly striking 
when writing the exact laws in MHD turbulence that 
stem from conservation of Et and He [IB] > &s mentioned 
before. As a result, the phenomenological formulation in 
Eq. (TTTJ)) seems the most plausible. 

If we now assume that Hc(k) scales as the total en- 
ergy, i.e. H c {k) ~ E T (k) ~ fc~ 3 / 2 , then Eq. (JT^ yields 
Hji(k) ~ k~ 2 . However, it should be noted that the spec- 
trum of cross-helicity in MHD turbulence may be steeper 
than IK, as, e.g., Hc(k) ~ fc _2 ,a spectral law which was 
obtained in two-point closures [45[ . In this case the spec- 
trum leads to H R (k) ~ k~ 5/2 . For 9 < k < 50, we find 
Ha(k) ~ A: -2 2 which does not allow us to decide between 
the two solutions. Moreover, note that the correlation co- 
efficient is very low in this computation, of the order of 
pc = Hq/Et ~ 1.5 x 10~ 3 , with a spectrum strongly 
fluctuating with wavenumber. 

These simple arguments indicate that there may be in 
fact dynamical constraints on the behavior of the flow 
that link the velocity and magnetic field statistics, both 
for the energy and the helicity, and that go beyond a 
strict derivation of scaling laws in turbulence based on 
the conservation of a sole invariant, e.g., the magnetic 
helicity or the total energy 



VI. DISCUSSION AND CONCLUSION 

In this paper we presented results from numerical sim- 
ulations of MHD turbulence in the absence of external 
forcing, using a Lagrangian averaged model to attain the 
largest possible Reynolds number at a given spatial reso- 
lution. At very high Reynolds number, the separation of 
scales between the energy-containing range and the dis- 
sipation range is huge and, with a power-law decrease of 
the distribution of energy among scales, for all practical 
purposes the large-scale magnetic field acts as a locally 
uniform field for the small scales; with one important dif- 
ference though: that there is no imposed direction to the 
flow and that it is in fact isotropic at large scales (no 
k = component). In other words, we can say that we 
have global isotropy and homogeneity. This is the case 
in the Solar Wind at very large scale: the locally uniform 
field in one of the spiral arms, changes direction with the 
spiral arm, and globally can be viewed as isotropic. How- 
ever, at small scale, the effect of the large-scale field is 
such that the dynamics is anisotropic, or more precisely, 
locally anisotropic. 

As a result of the global isotropy in our run, the to- 
tal energy spectrum develops an inertial range compati- 
ble with Iroshnikov-Kraichnan phenomenology, and other 
spectral quantities are consistent with this scaling. More- 
over, the small scales are seen to go through an Alfv- 
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enization process that can be attributed to a local mean 
magnetic field. 

However, magnetic helicity presents a wide inertial 
range with a spectrum steeper than what usual phe- 
nomenological arguments predict. What then is govern- 
ing the behavior of magnetic helicity? In the inverse 
cascade range, a dimensional Kolmogorov-like analysis 
leads to Hb(k) ~ k~ 2 (32|, and numerical data using the 
intcgro-differential MHD equations in the framework of 
the Eddy-Damped Quasi-Normal Markovian (EDQNM) 
closure does corroborate the existence of such a regime. 
One should note that the argument in [32j pre-supposes 
the independence of the cascades that may take place 
in MHD turbulence, that of total energy to small scales, 
that of cross-helicity likely to small scales as well, and 
that of magnetic helicity to large scales. However, the 
cquipartition relationships in Eqs. ([TU|) . (|T4")) . and (fT51) . 
link the different correlators that one can build on the 
physical variables in a much more intricate way. We do 
know already that there are other such complex relation- 
ships involving coupling of these quantities, namely the 
exact laws that stem from the conservation of Et, He 
and H b EElIzl. 

It is advocated in [35[ that, in the decay and the forced 
case, different dynamics may take place. It is true that 
for a steady-state inverse cascade to develop, separation 
of scales and forcing are in principle required, and the 
present simulation has initial conditions with most of the 
magnetic helicity concentrated at the largest scale in the 
box. At this point, it is therefore important to make a 
distinction between the behavior of magnetic helicity in 
its inverse cascade range, and its behavior in the direct 
cascade range of energy as observed in our simulation. 
In the direct cascade range, helicity seems to be trans- 
ported to smaller scales by the energy cascade, and its 
spectrum can be explained from phenomenological argu- 
ments that take into account Alfvenization of the flow, 
or equivalently, the dynamics of magnetic field induction 
leading to Alfvenization. 

These arguments are however dependent on the spec- 
trum of the kinetic and magnetic energy separately, on 
the spectrum of kinetic helicity H v , and on the spectrum 
of He- Indeed, the exact relationship for magnetic helic- 
ity derived in [47] is more complex in its structure than 
the laws for the other invariants; it involves the fields 
themselves, not structure functions (based on field dif- 
ferences). Moreover, the electromotive force v x b (asso- 
ciated with the magnetic induction) plays a central role 
(the law in [47| derives only from the induction equation). 
Furthermore, the third-order correlators of third-order 
tensors involving pscudo- vectors have a much more com- 
plex structure than for the full symmetric (non-helical) 
case, since it requires a priori four scaling functions to 
define the properties of the dynamics. Other examples 
where the presence of helicity leads to new scaling laws 
are known [46]. Thus, we can expect, in a similar fash- 



ion, that there could be different forms of the magnetic 
helicity spectrum in MHD depending on, e.g., the kinetic 
helicity or the intensity of the cross-helicity in the flow. 

These helicity effects open the door to a lack of uni- 
versality for spectra of invariants (see also [HI, [HI, 148[ ) , 
since different regimes may occur in different MHD flows. 
Note that these arguments differ from the critical balance 
regime for MHD turbulence advocated in [3l| and which 
does not seem to be observed in high Reynolds number 
MHD, at least in the globally isotropic and decaying case 
(see [15(, and Fig. 0|). In spite of the different possible 
scaling laws arising for magnetic helicity, it would be in- 
teresting to know if all satisfy the condition given by Eq. 
(fT8"| . We know of at least a few more examples satisfying 
this condition, e.g., the flow in Ref. [37|, and the flows 
studied in the pioneering work of Mtiller and Malapaka 

Finally, and as mentioned in the introduction, one 
other peculiarity of MHD turbulence is the fact that in- 
teractions are less spectrally local than in hydrodynamic 
turbulence. This has been known for a long time in 
its simplest version, that of Alfven waves. The large- 
scale magnetic field is responsible for the propagation of 
waves that put in equipartition the small-scale velocity 
and magnetic fields; this is the basis for the modification 
to a Kolmogorov energy spectrum as proposed in [28| . In 
fact, non-local effects have been quantified in MHD tur- 
bulence (see @, E3] and references therein) . Nonlocality 
may also break down the underlying hypothesis leading 
to Hb(k) ~ k~ 2 even for the inverse cascade range, and 
may intermingle small-scale velocity and magnetic exci- 
tations with the large-scale magnetic field. As stressed in 
[48| . nonlocality breaks the simple self-similarity known 
to occur in fluids at least at second order; it is also shown 
in [48| that the magnetic field at small scale may have 
a folded structure that is consistent with the exact law 
derived in [l6|. Thus, the exact MHD laws, and their 
associated balances presented here, which are derived in- 
dependently of locality assumptions, may therefore prove 
useful to reach a better understanding of MHD turbu- 
lence. 
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